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Solitary waves in the Nonlinear Dirac Equation with arbitrary nonlinearity 
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We consider the nonlinear Dirac equations (NLDE's) in 1+1 dimension with scalar-scalar self 

interaction ^i(**) fc+1 , as well as a vector-vector self interaction ^ I (*7 M *#7 M *)5 ( ' i:+1) . We 
find the exact analytic form for solitary waves for arbitrary k and find that they are a generalization 
of the exact solutions for the nonlinear Schrodinger equation (NLSE) and reduce to these solutions 
in a well defined nonrelativistic limit. We perform the nonrelativistic reduction and find the l/2m 
correction to the NLSE, valid when \u — m\ <C 2m, where lu is the frequency of the solitary wave in 
the rest frame. We discuss the stability and blowup of solitary waves assuming the modified NLSE 
is valid and find that they should be stable for k < 2. 

PACS numbers: PACS: ll.15.Kc, 03.70.+k, 0570.Ln.,11.10.-s 



I. INTRODUCTION 

Beyond the usual applications in field theory, the non- 
linear Dirac equation (NLDE) also emerges in various 
condensed matter applications. An important example 
being the Bose-Einstein condensate (BEC) in a honey- 
comb optical lattice in the long wavelength, mean field 
limit [ ]. The multi-component BEC order parameter 
has an exact spinor structure and serves as the bosonic 
analog to the relativistic electrons in graphcnc. 

Classical solutions of nonlinear field equations have a 
long history as a model of extended particles [_', 3]. The 
stability of such solutions in 3+1 dimensions was studied 
in detail by Derrick [ ] . He showed that the classical so- 
lutions of the self-interacting scalar theories (with both 
polynomial and non-polynomial interactions) were unsta- 
ble to scale transformations. However he was not able to 
make any conclusive statements about the spinor theo- 
ries. In 1970, Soler [3] proposed that the self-interacting 
4-Fermi theory was an interesting model for extended 
fermions. Later, Strauss and Vasquez [5] were able to 
study the stability of this model under dilatation and 
found the domain of stability for the Soler solutions. 
Solitary waves in the 1+1 dimensional nonlinear Dirac 
equation have been studied [ , ] in the past in case 
the nonlinearity parameter fe = 1, i.e. massive Gross- 
Neveu [8] (with N = 1, i.e. just one localized fermion) 
and massive Thirring [9] models). In those studies it 
was found that these equations have solitary wave solu- 
tions for both scalar-scalar (S-S) and vector-vector (V-V) 
interactions. The interaction between solitary waves of 
different initial charge was studied in detail for the S-S 
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case when fe = 1 in the work of Alvarez and Carreras 
[ ] by Lorentz boosting the static solutions and allow- 
ing them to scatter. Stability of the k = 1 problem was 
also studied by Bogolubsky [11], who found using a vari- 
ational method that preserved charge, that the frequen- 
cies u < l/v2 should be unstable. However, subsequent 
numerical work by Alvarez and Soler [12] showed that 
this result was incorrect (i.e. the solitary waves were nu- 
merically stable). Further analytic work on stability for 
the S-S model using the Shatah-Strauss formalism [13] 
by Blanchard et al. [14] turned out to give inconclusive 
results in that they could not prove that the solutions to 
the Dirac equation were minima of the variational energy 
functional. Thus the domain of stability of solutions to 
self interacting 4-Fermi theories is still an open question. 
In this paper we generalize the work of Lee, Kuo, and 
Gavrielides [ ] to arbitrary k and find exact solutions 
for all k. The paper is organized as follows: In Sec. II 
we find rest-frame solitary wave solutions of the form 
*(»,<) = e~ iut ip(x), for both the case of the S-S and 
V-V interactions. We calculate the rest frame frequency, 
w, and the energy, H , of a solitary wave of charge Q, as 
a function of the parameters k and g. We find the range 
of k and g values for which uj and H are in the range 
u 



0< 



H 



< m. In Sec. Ill we derive the nonrelativis- 



tic limit of the NLDE and find the leading term which 
is the nonlinear Schrodinger equation with corrections 
of the order of l/2m. Our derivation agrees with the 
heuristic result for k = 1 for modification of the NLSE 
found earlier by [15]. We find that the correction term 
has the same magnitude but opposite sign for the V-V 
as compared to the S-S case and find that the expansion 
is always valid whenever \u — m\ <C 2m. In the V-V 
case, the NLDE solutions are numerically quite close to 
those of the NLSE for all values of u. However for the 
S-S case, when we depart from the domain of validity of 
the non-relativistic reduction, the solitary wave solutions 
depart dramatically from the NLSE limit and become 
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double humped. We plot the crossover to this regime as 
a function of the nonlinearity parameter fc. In section IV 
we first discuss stability of solitary waves in the NLSE 
using an auxiliary Lagrangian for the static solutions. 
We find that the criteria for stability is < k < 2 and 
that identical results are obtained using stability against 
scale transformations (Derrick's theorem [4]). However 
the scale transformation argument leads to the conclu- 
sion that there should be unstable solitary waves in the 
NLDE for k > 1 which violates continuity argument to 
the nonrelativistic regime. It also led to contradictions 
with numerical experiments at k = 1. We then discuss 
the stability question in the modified NLSE (mNLSE) 
and show that it is essentially the same as for the NLSE. 
In section V we discuss how to obtain information about 
self-focusing in case k — 2 and k > 2 for both the NLSE 
and mNLSE assuming that the time dependent solitons 
are self-similar generalizations of the exact solution of the 
NLSE. We find that the correction terms in the mNLSE 
eventually dominate at late times during self-focusing 
and so the approximation breaks down during the late 
stages of self-focusing. 

We conclude with a summary of our main findings as 
well as a discussion about the possible future directions 
for settling issues of stability using various approaches 
including numerical methods. 



II. SOLITARY WAVE SOLUTIONS 

We are interested in solitary wave solution of the 
NLDE given by 

(i7^-TO)* + c/ 2 (M) fc * = 0, (2.1) 

for the scalar-scalar interaction and 

{il^df, - to)* + 5 2 7^*(*7 M *)(*7 At **7''*)5 (fc " 1) = , 

(2.2) 
for the vector-vector interaction. These equations can be 
derived in a standard fashion from the Lagrangian 



L = * (i7 M 9 M - to) * + Li . 
For scalar-scalar interactions, we have 



(2.3) 



L, 



-(**) 



fc+i 



fc + 1 
whereas for vector-vector interactions we have instead 



(2.4) 



fc + 1 



(2.5) 



Note that in the above equations, g 2 is the dimensional 
coupling constant, i.e. g 2 = G 2 m 1 ~ k , where G is dimen- 
sionless. The 7 matrices in 2 dimensions in our conven- 
tion satisfy 



{lmlv}+ = 2.9^; SV = 



1 
-1 



(2.6) 



We are looking for solitary wave solutions where the field 
* goes to zero at infinity. It is sufficient to go into the 
rest frame, since the theory is Lorentz invariant and the 
moving solution can be obtained by a Lorentz boost. In 
the rest frame we have that 



V(x 1 t)=e- tult ijj(x). 



(2.7) 



We are interested in bound state solutions that corre- 
spond to positive frequency in the rest frame less than 
the mass parameter to, i.e. < u> < to. For these bound 
state solutions one requires that the energy of the solitary 
wave H obeys < H < to. Choosing the representation 
7o = "3] 171 = <7i, where the <7j are the standard Pauli 
spin matrices, we obtain 

iosdt^ + <J x d x -$ - to* - V/* = 0, (2.8) 

where Vi = — Okf- . Defining the matrix, 



9* 



ip(x) 



m ubw cos ; 



(2.9) 



we obtain the following equations for u and v. For scalar- 
scalar interactions, we find: 



du 
dx 
dv 

dx 



+ (to + uj)v - g z (u A - v z fv = 



+ (to - uj)u - g z (u z - v z ) K u = 



(2.10) 



For the vector-vector case one has instead 

dv 
dx 
dv 
dx 



+ (m + u)v + g 2 (u 2 + v 2 ) k v = , 
+ (m - uj)u - g 2 (u 2 + v 2 ) k u = . 



(2.11) 



A first integral of these equations can be obtained using 
conservation of the energy-momentum tensor, 



T M „ = i^/j^d^ - g^L , 



d"T^ = ; 



(2.12) 
(2.13) 



which yields for stationary solutions 

Xio = constant , T\\ = constant . 
For all the cases we want to study we can write 

Tn = ui^tp - to^V + Li . (2.14) 

For solitary wave solutions vanishing at infinity the con- 
stant is zero and we get the useful first integral: 



Tn = tjip'tp — mtpip + Li = . 



(2.15) 



Multiplying the equation of motion for either the scalar- 
scalar or vector-vector interaction on the left by ip we 
have that: 

(k + l)Li = —uj^ip + rn^ip + Tpi~/idi?p . (2-16) 



We find from Eqs. (2.15) and (2.16) that 

ujkijj'i/j — mk'tpip + ^)i^\d\i\) = . (2-17) 

For the Hamiltonian density we have 

H = T 00 = i/ji"fidiip + mi/)ip — Lj 



hi + h 2 - h 3 



(2.18) 



Each of hi are positive definite. From Eq. (2.15) and 
(2.16) one derives that 



kLi = Tpijidiip 
which further implies that 

fr-3 = yhi . 



(2.19) 



(2.20) 



In particular, for fc = 1, we obtain % = rmpip. In terms 
of (R, 0) one has 



iftijidiip = %/j'x/j 



(19 
dx 



(2.21) 



This leads to the simple differential equation for for 
solitary waves 



cW 
dx 



= —LUk + rrik cos 29 , 



(2.22) 



where lu^ = kuj and m^ — km. The solution is 

9(x) = tan -1 (a tanh ^a:) , (2.23) 

where 



TOfc — ujk m — lu 



" = > ' = \l — ; — -■> Pk = J ml - lu 2 . 

m^ + LUk \ m + lu v" 

(2.24) 

In what follows it is often useful to rewrite everything in 

terms of a and /3. We have the relations: 

m + ui = — , to — lu = af3 , f3 = y to 2 — lu 2 . (2.25) 



A. Scalar-Scalar interaction 

First let us look at the S-S interaction. Using Eqs. (2.4) 
and (2.15) we obtain 



9 



. fc+i 



luR 2 - mR 2 cos 26» + -f— [R 2 cos 20) =0. (2.26) 

Thus 

'(k+ 1)(tocos20- to)' 



R 2 



g 2 (cos29) k + 1 



(2.27) 



We have 

d9 



$ 



dx Wfc + rrik cosh 2/3kX 



-wit + m k cos 20 , (2.28) 



so that 



cos 20 



?Tife + Wfc cosh2/3fea; to + w cosh 2/3^2; 



Wfe + TOfc cosh 2/3/ja; w + to cosh 2/3^2; 
One important expression is 



to cos 20 — w 



Using this we get 



$ 



fc 2 (w + to cosh 2f3]~x) 



(2.29) 
(2.30) 



2 a; + to cosh 2/3^2; 
to + w cosh2/3fea; 

Using the identities: 
1 + a tanh fi^x — 

1 — a tanh /3fc2; = 



(fc + lK 2 
g 2 k 2 (to + w cosh 2/3^x) 



(2.31) 



to cosh 2/3^2; + w 

m + lu 
w cosh 2/3^2: + to 

m + lu 



we obtain the alternative expression 



sech (3).x , 

sech flux , 

(2.32) 



2 1 + or tanh /3fcX 
R = 



1 — a 2 tanh /3fc2; 



(fc + l)^sech 2 /3 fc a; 



g 2 k 2 (m + u>)(l — a 2 tanh /3fcx) 

(2.33) 
The equation for lu in terms of g 2 is determined from the 
fact that the single solitary wave has charge Q, 



dx ijrij} 



dx R 2 {x) 



Thus the equation we need to solve for lu is 



1 



(* + !)# 



i/(fe) 



where 



/?fe [g 2 fc 2 (TO + cj) 



4[a 2 



1-a 2 ' 



(1 - y 2 )U k - l ){\ - aV)K*+i) 
For fc = 1, one obtains 

Ii [a 2 ] = 2 

Q = 

with the solution 

LU 



(2.34) 



(2.35) 



(2.36) 



(2.37) 



dx R 2 = 



4a _ 2/3 

(1 — a 2 )g 2 g 2 Lu 



v/l + QV/4' 



(2.38) 



in agreement with earlier results of [ ]. For fc = ^, we 
obtain 

» (l-y 2 )(l + a 2 y 2 ) 

J v (i-VF (2 - 39) 

(tanh - a — a) 



Iifa 2 ! = 



and 



q.MS" 2 , 



k 2 g 4 



For k — |, we obtain 



/ | [a 2 ] = -2if(a 2 ) 



4w (a 2 ) 

I -a 2 



(2.40) 



(2.41) 



where K (k) is the complete elliptic integral of the first 
kind. In general we can cast Ik [or] into the sum of two 
hypergeometric functions 2F1 ■ Letting y — x^ we have 
that 



h[a 2 



^x-ia + a^Kl-,)-^) 







From the definition 



rff 



^(l-t) 



:-&-! 



(l-a 2 x)i( fc+1 ) 

r(6)r(c-6) 



/o (l-te) a 

we find that 
.2 



r(c) 



2 Fi(a,6,c;^) , 
(2.43) 



h[a\k] 



(2.44) 



1 



111 
2'2 + fc' 

13 3 1 



which, when substituted into Eq. (2.35) is the equation 
we solve to obtain uj in terms of k, Q, m and g. Here 
B(x,k) denotes the Beta function. 

In order to see if the classical solution describes a 
bound state, one must calculate the value of the Hamil- 
tonian for this solution and show that it is less than m. 
The Hamiltonian density is given by (2.18), so that the 
energy of the solitary wave is given by 



dx % 



dx (h\ 



dx 



h ^ l -k 



h/l 2 

h 2 



h 3 ) 



= fli 1- 



#2, 



(2.45) 



where we have used Eq. (2.20). We find 

Mv) 



Hi = JdyR 2 (y) 

fa 



k(m + lu) 
x 



dy 



g 2 k 2 (m + lu) 



(2.46) 



B (h l+ l) 2 H l - 

H 2 = m / dy R 2 {y) cos 26 (y) 



1 1 3 
fe'2'2 



1 

Jk 

x B 



(fc + iM| " 

g 2 k 2 (m + u>) 

1 ^^f' 1 ' 



(2.47) 



2' k 



fc'2'2 



Without loss of generality, in the remaining part of this 
subsection we now put m — 1 so that < w < 1, in order 
to measure u, H in units of m. For k = 1 we find 



//, = 2(1 „ + ^ ) [(l + a 2 )tanh- 1 a -a], 



9 2 Q 



H, 



tanh a 



(2.48) 



Therefore we find that the energy of the solitary wave is 

4 



E„ 



9' 



where in a so i, uj so i = 



tanh' 1 a soh (2.49) 

We notice that the 



V!+Q 2 9 4 /4 

energy of the solitary wave with k = 1 does not depend 
on the width parameter /3. Simplifying we obtain for 
k = 1 and for all values of g 2 



H so i = dx if)ip — 



9 2 Q 



sinli (p 2 Q/2) < 1, (2.50) 



so that all the solutions are "bound states" . This agrees 
with the result of Lee et al. [6] . 



For k—\ one finds that 



Hi = 



H 2 = 



9(1 -u) 2 
16g 4 Q 

9(l + o;) 
2g i Q 



[{3a 4 + 2a 2 + 3) tanh 1 a - 3a(l + a 2 )] , 



[-a+ (1 + a 2 ) tanh 1 a] . 



(2.51) 



For Q = 1 and selected values of k we determine u> 
and -ffsoj and plot in Fig. 1 the allowed values for which 
H so i < 1. Note that the range of g values for the exis- 
tence of a bound state, as a function of A:, is bounded from 
below. The functional dependence of the lower bound 
9mim together with the corresponding solution uj(g m in), 
as a function of k, are depicted in Fig. 2. We note the 
rapid increase of g m in at large values of k. At k w 2, the 
upper bound of the solution u>(g m i n ) becomes lower than 
1, and we notice an inflection in g m in{k). Summarizing, 
we find that in the S-S case, bound states exist for all 
values of k and g > g m in- 



B. Vector- Vector interaction 

For the V-V interaction case, we obtain 

fc + i 



L, 



-Ji (\j/ 7AI ^y*)5( fc+1 ) = _JL_#2(/c+i) 



Eq. (2.15) now becomes 



(2.52) 



a 2 
)R 2 - mR 2 cos 26 + -JL_^R 2 ( k + 1 ) = . (2.53) 



Thus 



R 2 = 



(k + l)(mcos26 -lu) 



(2.54) 




FIG. 1: (Color online) NLDE bound states for the scalar- 
scalar interaction case: uj and H so i as a function of k and g 
for Q = l. 




FIG. 2: (Color online) Plot of the lower bound of the allowed 
range of g values in the scalar-scale interaction case, as a func- 
tion of k, together with the corresponding solutions, uj(gmin)- 
The solid lines are intended only as a guide to the eye. 



This can be rewritten in the following two forms: 
(k + l)(3 2 



It 



g 2 k 2 (uj + mcosh2(3i c x) 
(k + l)/3 2 sech 2 p k x 



(2.55) 



_g 2 k 2 (m + w)(l + a 2 tanh fax) _ 
The equation for uj can then be determined by using 



the charge defined in Eq. (2.34). This gives 



1 

fa 



g 2 k 2 (m + uj) 



i/(fe) 



Ik [a 2 



111 1 



where 

For k = 1, this gives 

4tan _1 a 



(2.56) 



(2.57) 



Q = 



dx It = 



9 1 



uj = mcos(g Q/2) 



(2.58) 
This imposes the restriction on the coupling constant, i.e. 
g 2 Q < it so that the spectrum is composed of positive- 
energy fermion states. On the other hand, for k = h this 
gives 



9(1+0;) _^ 2 

2j 4 



(2.59) 



The energy of the solitary wave is given by integrating 
the Hamiltonian density (2.18), and we obtain 



Hsoi = dx H = dx (hi + h 2 — h 3 ) 



^(l- k ) + H 2 



where 



Hi = 



fa 



k(m + uj) 



{k + l)Pl 
g 2 k 2 (m + uj) 



(2.60) 



(2.61) 



H 2 = 



1 
Tk 

x B 



(k + i)p 2 



113 1 

fc'2'2 + fc' 



g 2 k 2 (m + uj) 
1 1 



(2.62) 



U' J 



,111 1 

ill i 



lFl (rV 



Without any loss of generality, in the remaining part of 
this subsection we put Q = m = 1, i.e. we measure w, H 
in units of m so that < uj < 1 . For k — 1 , we find 



Hi = 

H 2 = 



2(1 \ ^ [a -(1- a 2 ) tan' 1 a], (2.63) 

g 2 

4a 



g 2 (l + a 2 ) ' 
For k = 1, we have an analytic solution: 

# soi = / dxtpip = — sin(.g 2 /2) < 1 



(2.64) 
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FIG. 3: (Color online) NLDE bound states for the vector- 
vector interaction case: ui and H so i as a function of k and g. 



since < g 2 < n, thereby showing the bound-state be- 
haviour even in the vector case. For k = |, one finds 

Hi = 9 ( 1 + ^ [(3a 4 - 2a 2 + 3) tan" 1 a - 3o(l - a 2 )] , 

g 2 = 9(1 + {j) [(l + a 2 )tan- 1 a- 



v " [( 



a(l — a 2 ), 
(1 + a 2 ) J 



#3 =2ifi 



In Fig. 3 we map out the allowed values of u and g 2 for 
various values of k. The allowed range of g values for the 
existence of a bound state, as a function of k, has both 
a lower and an upper bound, and the domain shrinks 
as k increases. Around fc=2.5, these bounds cross, and 
no bound states are possible for k > 2.5. The functional 
dependence of g m in and g m ax, together with the corre- 
sponding solutions u>(g m in) and uj{g ma x), as a function of 
k, are depicted in Fig. 4. As in the S-S case, co(g min ) be- 
comes less than 1 for k « 2, and we notice an inflection in 
gmin(k). However, we now find that io(g max ) approaches 
one in case k > 2. 



III. CONNECTION TO THE SOLUTIONS OF 

THE NLSE 

In this section we will perform the nonrelativistic re- 
duction of the NLDE to determine how it compares to 



FIG. 4: (Color online) Plot of the lower and upper bounds of 
the allowed range of g values in the vector-vector interaction 
case, as a function of k (top panel). In the lower panel we de- 
pict the k dependence of the corresponding solutions io(g m in) 
and uj(g max ), respectively. The solid lines are intended only 
as a guide to the eye. 



the NLSE. The NLDE can be written as 

ia 3 d t ^ + a x d x ^-m^-Vi^ = 0., (3.1) 

where Vj = — ^- = —g 2 (ip^) k . Next, we use Moore's 



(2.65) decoupling method [Hi] and write 



Vl[X] = l±^ Vl + X l^l Vl 



(3.2) 



We see that Vj [A = 1] = Vj . It has been shown that doing 
a perturbation theory in A is a valid way of obtaining the 
corrections to the nonrelativistic theory. Moore's decou- 
pling technique was used for the (relativistic) hydrogen 
atom using conventional Rayleigh-Schrodinger perturba- 
tion theory and computer algebra and it was shown that 
the perturbative solution converges to the correct solu- 
tion [ !]. It has been applied successfully to the rela- 
tivistic calculations on alkali atoms and represents one of 
the many relativistic perturbative schemes investigated 
by Kutzelnigg [ ]. We will show that this procedure 
leads to the heuristically derived nonrelativistic reduc- 
tion of the NLDE as discussed by Toyama et at for the 
case k = 1 [15]. 
We let 



*o(z) 



"0 

vo 



(3.3) 



be a solution of the theory when A = 0. For scalar-scalar 
interactions, we find: 



duo 
dx 
dv 

dx 



+ (to + uj)v = , 

+ (to - u)u - g 2 (uQU Q - VQv a ) k u = 



From Eq. (3.4) we obtain 

duo 
dx 



(to + ui)vq 



(3.4) 
(3.5) 



This leads to the following equation for uq: 



2to 



+ (Vi - co) (l + ^) u = , (3.6) 

where eo = lu — to. We notice that the expansion pa- 
rameter is £o/(2to). When \uj — m\/(2m) <C 1 is satisfied 
then we can be sure that the NLDE solutions go over to 
the NLSE solutions. However we will find that in the 
V-V case, the reduction numerically appears valid over a 
wider range. The relevant Schrodinger-like equation is: 



(uo)xx 





2m 


- + ViUq = Eu 


where 






Vi 


= trfl + M, 


E = e (l 



2to 



(3.7) 



(3.8) 



For consistency we need to expand Vj to first order in 
l/2m. For the scalar scalar case, we have 

Vi = -g 2 (u* u a -VQV ) k 



K«o) fe 



4m 2 



(«o«o) (uo)%(u )a 



(3.9) 



The resulting modified nonlinear Schrodinger equation 
(mNLSE) can be derived from the Lagrangian: 



L = iip*d t ip 

P. 2 



1 

2ro 



<■■>..■ (l + J^VO* 



fe+i 



k + 1' 
and the Hamiltonian is given by 



(3.10) 



H = 



dx 

2ro 



9 



2to 



9 



1 



(V»V) 



fc+1 



(3.11) 
where g 2 = g 2 [l + eo/(2m)]. 

In the case of the V-V interaction, the nonrelativistic 
reduction of the NLDE is similar to the previous case 
with the difference that 



Vj v 



-g 2 (uou 



«o u oJ 



K u o) fc + -r-^K«o) MxMx 



4to 2 



(3.12) 



The resulting modified nonlinear Schrodinger equation 
(mNLSE) can be derived from the Lagrangian: 



L = itp*d t ip 
9 2 



1 
2m 



$1$* 



1 



9 
2to 



(^» fe 



WW 



fc+1 



fc + 1 

and the Hamiltonian is given by 



(3.13) 



H 



dx 
2ro 



^x 1 



9 

2to 



(^V) fe 



1 



(V-V) 



fc+1 



(3.14) 
where <p = g 2 [l + £ /(2to)]. 

Thus we see that the resulting theory in the large 1/2to 
limit (as well as when \ui — m\ <C 2m), in both S-S and 
V-V cases reduces to the modified NLSE equation. The 
first correction has the same magnitude but opposite sign 
for the two cases. 



A. Comparison with the exact solution of the 

NLSE and mNLSE 



Here we want to compare the NLDE with the exact 
solution of the of the NLSE as well as mNLSE for arbi- 
trary fc. We will give numerical comparison both when 
the criterion |w — m\ ^C 2to is satisfied and for general w. 
We will find that the V-V NLDE case has solutions that 
track those of the NLSE for a broader range of lo. 

First let us obtain solutions to the NLSE for arbi- 
trary fc. The NLSE is defined by the Lagrangian 



L = 



dx (rA - r^) - h . 



where for the S-S interaction 

I 



H 



dx 



-\7ip*\7tp - g- 



(P1>) 



fc+i 



2m k 

This leads to the equation of motion 

. dip 1 ( dijj s 



yw>*V)> = o 



(3.15) 



(3.16) 



(3.17) 



dt 2to \ dx 
If we make the ansatz 
ip(x,t) = r(y) exp[i(mvy — ut + 8)] , y = x — vt, (3.18) 
then it is easy to show that r(y) satisfies the equation 

r"(y)~nr(y)+g 2 r 2k+1 (y) = 0, (3.19) 

Equation (3.19) has an exact 



where $7 
solution 



provided 



r(y) ■■ 



mi- 
2 







D 



Asech^ k [D(y + y )], 
2k _(k + l)D 2 



2mk 2 



A 



2mg 2 k 2 



(3.20) 



(3.21) 



The mass density in the rest frame (v = 0) is given by 



ip*tp 



(k + l)D 2 
2mg 2 k 2 



l/ft 



sech 2/k [D(x + x )} . (3.22) 



Let us now obtain the solutions of the mNLSE. We first 
notice that to the first order in l/2m, the static mNLSE 
equation in both S-S and V-V cases is given by 

-(u ) xx + (m 2 -u> 2 )u -(m+uj)g 2 [uQUo} k u Q = 0, (3.23) 

which has the exact solution 

u (x) = Asech 1/k [j3 k (x + x )] , (3.24) 

with 

2k _ (k+l)Pl 



A 



(m + oj)g 2 k 2 



Hence for mNLSE, the mass density in the rest frame 
(v = 0) is given by 



p = ip*ip 



(rn + uj)g 2 k 2 



sech 2/fe [/3 fc (x + a;o)] 



(3.26) 

We will now compare the NLSE and mNLSE solutions 
with the solutions of the NLDE. In making these com- 
parisons we will in all cases compare the solutions for the 
charge density (which is the mass density for the NLSE 
and mNLSE). 



0.20 



0.15 . 



(3.25) & 




FIG. 5: (Color online) Comparison of the NLSE and NLDE 
solutions in the case of scalar-scalar interactions for k — 1, 
and w/m = 0.9 (top panel) and ui/m — 0.3 (bottom panel), 
respectively. 



B. Scalar-Scalar interaction 



One can rewrite the charge density p = R 2 , Eq. (2.33) 
in the following form which isolates the previous solution 
to the NLSE. 



1) 



/3 fc 2 (fc 
g 2 k 2 (m + lo) 

I 



l/fc 



sech 2 ' /3feX f(a,(3,x), 



f(a,P,x) = 



a 2 tanh /3kX 



(f-a 2 tanh 2 ^ fe a;)( 1 + 1 /fc) 



(3.27) 



If we compare NLSE and S-S case, we find that p(x = 
0) is same in both cases only if we can identify D with 
(3 k . We also have that f(a,/3,x = 0) = 1, so that with 
this identification, the charge and mass densities have the 
same value as a function of k for the NLSE and NLDE. 

On the other hand, p(x — 0) is strictly identical for 
S-S and mNLSE cases and no identification needs to be 
made. 

We have seen that the nonrelativistic limit is obtained 
when \u — m\/2m <^ 1. In Fig. 5, we compare the so- 
lutions to the NLSE and NLDE when w/m = 0.9 (top 
panel) and Lu/m = 0.3 (bottom panel), for k = 1. In the 
latter case, we notice that the solution to the NLDE is 
double humped. For any u < uj c (k) for which the solu- 
tion becomes double humped in the NLDE is shown in 
Fig. 6. 



M 

s 




FIG. 6: (Color online) Critical value, U) c (k), for any ui < u) c (k) 
the solution of the NLDE equation becomes double humped 
in the case of scalar-scalar interactions. 
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FIG. 7: (Color online) Comparison of the NLSE and NLDE 
solutions in the case of vector- vector interactions for k = 1 
and cj/m = 0.01. 



C. Vector- Vector interaction 

Now we rewrite the solution found for the charge den- 
sity p = R 2 , Eq. (2.55) in the following form: 



/\ 2 (fc + l) 



l/fe 



sech 2/fc f3 k x f(a,0,x) 



_g 2 k 2 (m + u>) 
/(a, /3, x) = (1 + a 2 tanh 2 f3 k x)- 1/k 



(3.28) 



We have seen that the nonrelativistic limit is obtained 
when \lu — m\/2m <C 1. For the V-V case, the modifica- 
tion of the NLSE result is small even at very small uj/m 
and, unlike in the case of S-S interactions, the NLDE 
solution never becomes double humped. In Fig. 7 we 
compare the solutions of the NLSE and NLDE when 
uj/m = 0.01 and k = 1. The main difference compared to 
the S-S case is that the convergence to the nonrelativistic 
limit as tu/rn — > 1, occurs from above in the vector case 
instead of from below as in the scalar case (see Fig. 7 
and the top panel of Fig. 5). Again notice that p(x = 0) 
is identical in NLSE and V-V case only if we identify D 
with /3fc. 

On the other hand, p(x — 0) is strictly identical for 
mNLSE and V-V case and no identification needs to be 
made. 



IV. STABILITY OF STATIC SOLUTIONS 

The stability of the solitary waves of the NLSE have 
been studied for a long time. A recent discussion of this 
is found in [18]. In this section we will first show that 
an analysis of the solutions of the NLSE equation using 
the slope criterion ( ^J ' < for stability) where M 
is the mass of the Solitary wave and w the frequency 
gives the same result (0 < k < 2) as an analysis based 



on whether a scale transformations raises or lowers the 
energy of the solitary wave. The latter criterion is similar 
to the arguments first used by Derrick [4] in his study of 
the relativistic scalar field theories. We will then use a 
similar scaling argument first made by Bogolubsky [11] 
for the NLDE equation to obtain a criterion for stability. 
We will find that the results of this approach do not agree 
with a smooth continuation of the result for the NLSE. 
We will discuss the most likely reason for the failure of 
this method when applied to the NLDE. Finally we will 
look at the stability in the mNLSE which contains the 
first relativistic correction to the NLSE and show that it 
gives essentially the same criterion as that found for the 
NLSE, i.e. when < k < 2 we expect the solutions to be 
stable. 

Most studies of the stability of static solutions of the 
NLSE rely on the existence of a variational principle 



S£ = 5{H - ujM) = 



(4.1) 



from which the ordinary differential equation for the solu- 
tion u(x, oj) can be derived. Here the NLSE Hamiltonian 
is 



H = dx 



^WW-^M) 



fc+1 



and the mass is given by 



M 



dxij>*ip 



> (4.2) 



(4.3) 



This variational principle is quite similar to the one 
used to study the stability in the generalized KdV sys- 
tems [19-22]. There one derives the solitary wave equa- 
tion from 



Se = 5(H - cP) = , 



(4.4) 



where c is the velocity of the solitary wave while the 
generalized KdV equation is 



u t + u 2 u x 



+a[2u p u xxx + Apu p 1 u x u xx + p(p - l)u p 2 (u x ) 3 ] = . 

(4.5) 



This can be derived from the Hamiltonian 



H 



dx 



au p (u x )' 



1(1-1) 
and the corresponding momentum P is given by 

dx-u 2 (x,t) . 



(4.6) 



(4.7) 



Stable solitary waves of the form ip(x, t) = u(x, uj)e~ lut 
need to be local energy minimizers of the functional (4.1). 
Based on linearized perturbation theory and using this 
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variational principle Vakhitov and Kolokolov [21] showed 
that a necessary criterion for stability is that 



dM{uj) 

did 



<0. 



(4.8) 



This criteria is the analogue of the result found for the 
generalized KdV equations by Karpman [19] and Dey 
and Khare [20] who obtained that stable solitary waves 
for that system of equations required 



dP{c) 

dc 



>0. 



(4.9) 



The exact solution for u(x) of NLSE for arbitrary k is 
given in (3.20). Using that solution, one finds that the 
mass has the following dependence on ui: 



M = C(- w )( 2 - fe )/< 2fe ) 



(4.10) 



and the necessary criterion for stability is 

k < 2. (4.11) 

Another approach to stability, which leads to the same 
result as (4.11), is based on whether a scale transforma- 
tion which keeps the mass M invariant, raises or lowers 
the energy of a solitary wave. For the NLSE with Hamil- 
tonian 



H 



dx 
H\ — H 2 



d x ii)*d x ii} 



!J 



k + 1 



(V>V) 



fc+i 



(4.12) 



both Hi and H 2 are positive definite. A static solitary 
wave solution can be written as 



ip{x,i) = r{x)e 



(4.13) 



The exact solution has the property that it minimizes 
the Hamiltonian subject to the constraint of fixed mass 
as a function of a stretching factor (3. This can be seen 
by studying a variational approach as done in [22] or by 
directly studying the effect of a scale transformation that 
respects conservation of mass. 

In the latter approach, which generalizes the method 
used by Derrick [ ], we let 



and consider 



f3x 



ipp{x) = P 2 r{j3x)e 



this leaves 



M 



dx Tp*tp = I dx l/jalp/3 



(4.14) 



(4.15) 



(4.16) 



unchanged. One defines Hp as the value of H for the 
stretched solution ipp. One then finds that 



dH B 



d& 



/3=1 



0, 



(4.17) 



is consistent with the equations of motion, and the stable 
solutions satisfy 



d 2 H p 

d/3 2 



>0 



(4.18) 



If we write H in terms of the two positive definite pieces 
H\, H 2 , then 



Hp = f3 2 Hi - f3 k H 2 



We find: 



We obtain 



dHp 
d/3 



2[3Hi - {k)p k -"H 



(4.19) 



(4.20) 



^\p =1 = ^ Hl = ^H 2 . (4.21) 

This result is consistent with the equation of motion. The 
second derivative is given by 



d 2 Hp 
d/3 2 



2Hi - k{k - \)^~ 2 H 2 



(4.22) 



which when evaluated at the stationary point yields 



d 2 H p 

d/3 2 



= 2(2 - k)Hi 



(4.23) 



This result indicates that solutions are unstable to 
changes in the width (compatible with the conserved 
mass) when k > 2. The case k = 2 is the marginal 
case where it is known that blowup occurs at a critical 
mass (see for example Ref. 22). The result found above 
for the NLSE has also been found by various other meth- 
ods such as linear stability analysis and using strict in- 
equalities. Numerical simulations have been done for the 
critical case k = 2 showing that blowup (self-focusing) 
occurs when the mass M > 2.72 [23]. For k > 2 a vari- 
ety of analytic and numerical methods have been used to 
study the nature of the blowup at finite time [24] . 

Let us now apply this scaling argument, as was done by 
Bogolubsky [ ]), to the 1+1 dimensional NLDE. Again 
we will assume that the exact solution minimizes Hp 
when (3 — 1 with the constraint that the charge is kept 
fixed. (The validity of this assumption will be challenged 
below. All that is known is that Hp is a stationary point 
at the solution.) 

Our exact solution is of the form 

Because we want to keep the charge fixed, we consider 
the following stretched solution: 

^)=(:)=Mwx)(-w) e -«. (4 . 25) 
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The value of the Hamiltonian 

H = I dx -tjji'yiditp + rwipip — — 
= Hi + H 2 - H 3 , 
for the stretched solution is 

Hp = pH x + H 2 - p k H 3 



■(H>) 



fe+1 



(4.26) 



(4.27) 



where again H t are all positive definite. The first deriva- 
tive is 



9 ^=H l - k^H 3 
dp 

At the minimum, setting f3 = 1, we find in general 

1 



(4.28) 



H 3 = -H, 



(4.29) 



which is consistent with the equation of motion result we 
obtained earlier, see Eq. (2.20). We see that for fe = 1 
the energy is given by just H 2 . The second derivative 
yields: 



d 2 H, 
dp 2 



P 



-k{k-l)^- 2 H 3 . 



(4.30) 



From this we see that if k > 1, this analysis (if correct) 
would suggest that solitary waves are unstable to small 
changes in the width. For < k < 1 the solitary waves 
are stable to this type of perturbation. This argument 
does not depend on L/ as long as L/ is positive definite. 
The same result is valid for both scalar and vector type 
interactions. 

For k = 1, this argument does not give any insight 
into whether the solutions are stable. However, it is 
known that the solitary waves discussed here for k = 1, 
do appear to be stable numerically. Further, when they 
are scattered in numerical experiments, they interchange 
charge and energy, and sometimes show bound state pro- 
duction. Detailed numerical simulations have been per- 
formed by Alvarez and Carreras [10]. These results con- 
tradict the work of Bogolubsky [11] who studied changes 
in the frequency u> while keeping the charge fixed. There 
a similar analysis gave a maximum for the Hamiltonian 
when u < l/\/2, even though, as remarked above, nu- 
merical studies show that solitary waves in that frequency 
range are in fact stable. 

We have already shown above that the solutions of the 
NLDE reduce to those of the NLSE in the nonrelativistic 
limit. Assuming continuity arguments apply, one would 
expect that there would be at least a range of values of ui 
for which the solutions to the NLDE are stable for k < 2. 

So one needs to understand the reason for this apparent 
discrepancy. The main reason for assuming instability 
when the second derivative of H(/3) is positive, is that the 
stable solutions to the Dirac equation are at least relative 
minima of the effective action. However, the study by 
Blanchard et al. [14] to find an analytic criterion for 



stability in the 1+1 dimensional NLDE using the Shatah- 
Strauss formalism found that bound states were not local 
minima on the manifold of constant charge. This result 
is quite different from what happens in the NLSE where 
the bound states are local minima on the manifold of 
constant mass. So one cannot assume that the sign found 
in Eq. (4.30) yields information about the stability of the 
solution. On the other hand we can assume by continuity 
that there is a region where the analysis of stability in 
the mNLSE will give us information about stability at 
least in the regime where the expansion parameter e/2m 
is small. For the mNLSE we can use the scaling argument 
or the auxiliary variational approach to discuss stability. 
It is interesting that Derrick [4] in his seminal paper was 
unable to find a suitable method for discussing stability 
for self-interacting spinor theories. 

For the mNLSE the Hamiltonian for the S-S interac- 
tions is given by 



H 



dx 
2m 



Tpl^x 



1 



2m 



(W)' 



9 
fc+1 



W$) 



fc+1 



It is well known that using stability with respect to scale 
transformation to understand domains of stability applies 
to this type of Hamiltonian. This Hamiltonian is a sum 
of two positive and one negative term i.e. 

H = H 1 +H 2 - H 3 . (4.32) 

For the V-V case, the Hamiltonian is instead 

H = H 1 -H 2 -H 3 . (4.33) 

We also know that H 2 is of order g 2 /2m and is presumed 
small. If we again make a scale transformation on the 
solution which preserves the mass M — J ip*ipdx, 



we obtain 



4> f3 = p 1 ' 2 ^) , 



H = f3 2 H 1 ± f3 2+k H 2 - f3 k H 3 



(4.34) 



(4.35) 



Here the upper (lower) sign corresponds to the S-S (V-V) 
case. The first derivative is: 



dH 



2/3H 1 ± (2 + k)(3 k+1 H 2 - kp k ~ x H 3 . (4.36) 



Setting the derivative to zero at = 1 gives the equation 
consistent with the equations of motion: 

kH 3 = 2Hi ± (2 + k)H 2 . (4.37) 

The second derivative at /3 = 1 can now be written as 

d 2 H 



d(3 2 



= (4 - 2k) Hi ± 2(2 + k)H 2 



(4.38) 



This will be positive for k < 2 and the addition of a small 
H 2 should extend the stability of the solutions beyond 
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k = 2 in the S-S case. However, in the V-V case there is 
a somewhat lower region of stability. At k = 2, as we shall 
see below, the usual NLSE solitary waves blow up once 
the mass exceeds a critical value. For k = 1, numerical 
experiments for the time evolution of an initial wave of 
the form 

ip(x, t = 0) = ^/J3j2sech(l3y)e i{mvx ~^ mv2t - eot) (4.39) 

at t = relaxed to an exact solitary wave solution of 
the mNLSE that was not very different than the NLSE 
solution [15]. This result supports the conclusion that 
the solitary waves of the mNLSE are stable for k=l. 



V. SELF-SIMILAR ANALYSIS OF BLOWUP 
AND CRITICAL MASS FOR THE NLSE AND 

THE MNLSE 

To study in a "mean field" approximation blowup and 
critical mass, we look for self-similar solutions of the 
form: 

ip(x,t) =A(t)f(j3y)expi [mvy + A(t)y 2 + ivt] . (5.1) 

Here A(t),A(t) and f3(t) are arbitrary functions of time 
alone, and y — x—vt. What we have in mind is to start at 
t = with the exact solution of the form A sech 1 ' (Dy) 
and assume that this solution just changes during the 
time evolution in amplitude and width conserving mass. 
With this assumption one can derive the dynamical equa- 
tions for A and D from the action principle. The action 
for the NLSE is given by 



dtL 



where L is given by 



i = - 1 d d x{r^ t - vavo - h , 



with 



H = dx 



Vx^x 2^*^) 



*J,\k + l 



2m 



k + 1 



(5.2) 



(5.3) 



(5.4) 



The NLSE follows from the Hamilton's principle of least 
action: 



ST 8T 



(5.5) 



The NLSE has three conservation laws: mass, momen- 
tum and energy which can be derived from Noethcr's 
theorem in the usual fashion. The conservation of mass 

f A 2 f°° 

M = / ip*tpdx = —d, C 1= f(z)dz , (5.6) 

J P J — oo 



allows one to rewrite A(t) in terms of the conserved mass 
and the width parameter /3 and a constant C\ whose 
value depends on f{z). Thus, 



l2 M(3 



A 1 = 



C\ 



For f(z) = scch 7 (z), one obtains 

/^r( 7 ) 



Ci = 



r(7 + J) 



(5.7) 



(5.8) 



First consider the kinetic energy (KE) term in the La- 
grangian Density 






mv — Ay + 2vAy — ui 



(5.9) 



Integrating over space and scaling out j3, we obtain 

(5.10) 



KE/M = rav 2 - u> - AG 2 -?- 



where G — 4 and 



P 



Co 



z 2 f 2 (z)dz 



(5.11) 



■4 7 1 4^3(7,7,7,27; 7 + 1,7 + 1,7+1;-!) 



Next consider 



We obtain 
H /M 

where 

C 3 - 



H 



dx 



1 dip* dip 
2m dx dx 



(5.12) 



^ + ^^+4A 2 ^ (5 13) 
2 + C 1 2mG 2+4A C 1 2m' (5 ' 13) 



°°ifm dz -^ T{ " +1) 



2r( 7 + l) 



Finally for the interaction term: 



Hr = - 



fc + 1 



we obtain 

Hl/M = 

where 

c\ = / f 2k+2 Hz)dz 



dx(ijj*ip) 



Ca ( M 



*„i,\k+l 



(fc + l)Ci \C 1 G / 

v^r[(fc + i) 7 ] 



r[(fc + i) 7 - 



(5.14) 



(5.15) 



(5.16) 



(5.17) 
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Putting this together we get the following "effective La- 
grangian" for the time dependent functions G, A: 



L = 



-lo-AG 



2C2 C3 



1 



d d 2mG 2 



(5.18) 



4A : 



2 C 2 G 



9 



C*( M \k 



d 2m (k + 1) 
Lagrange's equation for A yields 



O4 / M \ 

Y)d\dG) 



A 



2mG 
AG 



(5.19) 



The first integral of the second order differential equa- 
tion resulting from the Lagrange's equation for G can 
be obtained by setting the conserved Hamiltonian to a 
constant E. One then has 



E 



C, 



1 



4A : 



2 C 2 G 2 



Ca 1 M 



Af*_V 

\\dGJ 



d 2mG 2 d 2m (k + 1) d 

(5.20) 
Using Eq. (5.19) we obtain the first order differential 
equation for G: 



E = 



C 2 2mG 2 C 3 1 g 2 C 4 / M 

d 4 + d2mG 2 ~ (k + l)d\dG' 

(5.21) 
We notice that at the critical value of k = 2, that the last 
two terms both go like 1/G 2 . Self- focusing occurs when 
the width can go to zero. Since G 2 needs to be positive, 
this means that at k — 2, the mass has to be greater than 
or equal to M* for G to be able to go to zero. Here 



M* 



or 



f 2m~igM k 



/ M " 



l 3C 2 C 3 

c 4 



dj_ 

C 4 2m ' 



1^ 



2.7207. 



(5.22) 



(5.23) 



provided we use the exact solution for k — 2, namely 
/ = sech 1 "(z) (which is a zero-energy solution). This 
agrees well with numerical estimates of the critical mass 
[23] and is slightly lower than the variational estimate 
obtained earlier by Cooper et al. [ ] using post-Gaussian 



trial wave functions. In the supercritical case we have 
that 



C 2 2mG 2 
d 4 



C A / M 



(fc + 1) 



»£(<&)■ <-) 



Thus G approaches zero in a finite time in this self-similar 
approximation with critical index: 



G^(t- t c ) 2 /( fc + 2 ) 



(5.25) 



This "mean-field" result was obtained earlier in [22, 25]. 
Now we would like to see how this argument is modified 
when we add the J^- corrections coming from the non- 
relativistic reduction of the NLDE. We now have: 



L 



dx(ip*ip t - ip* t ip) - H . 



(5.26) 



where for the mNLSE, the Hamiltonian is given by 



H = I dx 
9 



1 

2m 

2 



lm 



(5.27) 



k + 1 



W» 



fc+i 



Here upper (lower) sign corresponds to the S-S (V-V) 
case. Now we get one more term in the energy conser- 
vation equation. Also Lagrange's equation for A gets 
modified. The new term is 



5H/M = ± 
= ± 
x 
where 
Ex-. 

and 



1 



M Am 2 



dx^x (^» fe 
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M 



Am 2 \d 
Ci 



/\2 f (2fe+2) 



c 2 

d' 



i 2 v 2 G~ k 



Eh 

C x 



(5.28) 



2/^i2 — k 



AA 2 G 



(f'Tf 



(z)dz 



/^ 7 2 r[(fc + 2)7] 



2r[(fc + 2) 7 +|] 



(5.29) 



22(k+i)~t-i 
E-j = I z 2 /( 2fe + 2 ) ( z )dz = T — — . 4^3(^7 + 7, A;7 + 7,fe7 + 7,2A;7 + 27;fe7 + 7 + l,A;7 + 7 + l,fc7 + 7 + 1;-1) 



(fc + l) 3 7 3 



(5.30) 



Lagrangian's equation for A now yields 



A = 2m 



G_ 

AG 



l^fflV 



2m \dj C; 



Conservation of energy in the comoving frame (v = 0) 



(5.31) 



now leads to 

C 3 1 



E 



d 2mG 2 



4A 



d 2m 



9 2 C 4 



± 



9 



Am 2 



M 



El 



(k 

C~i 



l)Ci 



G -(fc+2) + ^ 4A 2 G 2-fc 



M 

, (5.32) 



or 



E=- 



C, 



1 



d 2mG 2 



1± 



.g 2 Ex f M 



2m C 3 VCiG 



2m /,2 C '2 

~T G c7 

(fc + 1) Ci 



1± 



<y 



M 

cTg 



E'2 

2m~C~ 2 

k 



M 
C\~G 



From this expression we again see that k = 2 is the crit- 
ical value. If the initial value of G is large enough so 
we can ignore the g 2 /2m corrections then in order for 
G 2 > 0, so that the width can decrease, one needs that 



f7 hngM* > 



I3C 2 C 3 



C 4 



(5.33) 



When G gets very small then the g 2 /2m corrections get 
large and our expansion breaks down. Blowup then needs 
to be studied using the full NLDE. We intend to do nu- 
merical studies of blowup in the NLDE in the near future. 
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a variety of ways. By continuity we expect that at least 
in the regime where the solutions of the NLDE are small 
perturbations of those of the NLSE, the solutions we have 
found will be stable for k < 2. We discussed the case 
k = 2 for the mNLSE approximation in detail as well as 
blowup for k > 2 using a self-similar ansatz. 

Before ending wc point out some of the possible open 
questions. 

1. Is there a connection between instability and the 
double hump behavior? 

2. In the V-V case we notice from Fig. 4 that while for 
k < 2, u(g m i n ) > Lo(g max ), for k > 2, the opposite 
is true. Is this somehow related to the fact that the 
NLDE V-V bound states are stable (unstable) for 
k < (>)2? Further, the dip in the value of ui(g m i n ) 
precisely occurs around k = 2 in both the S-S and 
the V-V cases. Is that just a coincidence or is it 
related to the instability for k > 2? 

3. For k = 1, it is known that the bound states of N 
localized fermions are stable in both the S-S and 
V-V cases. It would be interesting to examine if 
this continues to be true for arbitrary positive k. 

We hope to address some of these questions in the near 
future. Also we intend to do numerical simulations of 
collisions to see how energy and charge are exchanged, 
and also study blowup to understand whether there is 
much difference between self-focusing in the NLDE and 
the NLSE. 



In this paper we have found new solutions to the NLDE 
with arbitrary nonlinearity parameter k in the case of 
both the S-S and V-V interactions. The solutions for the 
S-S interactions have the property that for 10 > u> c (k) the 
shape of the solitary wave is similar to a sech 7 (x) pro- 
file, whereas for lo < Lo c {k), the shape is double humped. 
In the V-V case, the shape of the profile is always of 
the form sech 7 (x). We discussed the nonrelativistic re- 
duction of the NLDE and obtained a modified NLSE 
(mNLSE) whose stability properties could be studied in 
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